data <- read.csv(file="data1.csv",header=FALSE)
rows <- dim(data)[1]
columns <- dim(data)[2]
column_titles <- as.vector(t(t(data[1,2:columns])))
row_titles <- as.vector(t(t(data[1:rows,1])))
colors <- rainbow(columns)
par(bg="gray90", mfrow = c(rows-1, 1))
for(i in 2:rows){
 data_tmp <- as.matrix(data[i,2:columns])
 r <- barplot(data_tmp, names.arg = column_titles, col=colors, sub=row_titles[i])
 if(i == 2){
 title(main="Coreactants choice, A=1E5, Ac=1E7, O=0.01, JC")
 }
}
